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We investigate the properties of eigenstates and local density of states (LDOS) for a periodic 
2D rippled billiard, focusing on their quantum-classical correspondence in energy representation. 
To construct the classical counterparts of LDOS and the structure of eigenstates (SES), the effects 
of the boundary are first incorporated (via a canonical transformation) into an effective potential, 
rendering the one-particle motion in the 2D rippled billiard equivalent to that of two-interacting 
particles in ID geometry. We show that classical counterparts of SES and LDOS in the case of 
strong chaotic motion reveal quite a good correspondence with the quantum quantities. We also 
show that the main features of the SES and LDOS can be explained in terms of the underlying 
classical dynamics, in particular of certain periodic orbits. On the other hand, statistical properties 
of eigenstates and LDOS turn out to be different from those prescribed by random matrix theory. 
We discuss the quantum effects responsible for the non-ergodic character of the eigenstates and 
individual LDOS that seem to be generic for this type of billiards with a large number of transverse 
channels. 

PACS: 05.45+b, 03.20. 



I. INTRODUCTION 



The goal of this study is to deepen our understanding of the quantum-classical correspondence for chaotic systems 
by analyzing the properties of eigenstates and local density of states (LDOS) in the energy representation. This novel 
approach has already been successfully applied to models of two interacting particles jl|-|3| and to a three orbital 
schematic shell model [Q. In this paper we extend this approach to chaotic periodic billiards by incorporating the 
effects of the boundaries into the Hamiltonian operator and treating the degrees of freedom as independent "particles" . 
Some results of this work were advanced in j^] . 

The outline of the paper is as follows. Section II presents some basic features of the classical dynamics of the rippled 
billiard, together with the necessary details of its quantum description. In Section III the classical representation of 
the model is given. In Section IV we discuss the meaning of the structure of eigenstates (SES) and LDOS, and in 
Sections V and VI we define their classical counterparts. In Section VII we compare the quantum and classical SES 
and LDOS. Section VIII pertains to the individual properties of eigenstates and LDOS. In Section IX we make some 
clncluding remarks. 



II. DESCRIPTION OF THE MODEL 



The system that we shall use to explore the quantum-classical correspondence is a 2D billiard with periodic boundary 
conditions along the longitudinal coordinate x. The top profile is given by the function y = d + a^(x) with £(x + 2tt) = 
where a is the amplitude of the ripple and d is the width of the billiard when a is zero. The bottom boundary 
is flat (see Fig. 1). 

The first studies of the classical dynamics of this system were performed in Ref. Q in connection with the beam- 
beam interaction problem (see also 0). More recently, the finite length version of this system was analyzed in Ref. ||] 
as a model of a mesoscopic electron waveguide, where a transport signature of chaos in the behavior of resistivity was 
established. The same model has been also used |J to study certain quantum transport properties through ballistic 
cavities. On the other hand, the analysis of the band-energy spectrum for an infinitely long rippled channel |lC| ] in 
the cases of mixed and global classical chaos, provides insight into the universal features of electronic band structures 
of real crystals JO] . In Ref. |12| a similar periodic billiard was considered in connection with the problem of chaotic 



diffusion in the systems with band spectra, see also Refs. Jl3 14|. 

An attractive feature of this rippled channel is that its classical phase space undergoes the generic transition to 
global chaos as the ripple amplitude a increases. Hence results obtained for this particular system are expected to be 
applicable to a large class of systems, namely non-degenerate, non-integrable Hamiltonians. 

In the quantum description the model is given by the Hamiltonian 

H = ^-(P 2 + P 2 y ) =-^~{dl + d 2 y ) (1) 

for the wave function ^{x,y) obeying the boundary conditions, ^>(x,y) = 0, at y — and d + a£(x). There exist 
various numerical methods that can be used to obtain the eigenvalues and eigenfunctions of non-integrable billiards, 
such as the transfer matrix approach the scaling method p6| , and the scattering approach Jjl|. Here we shall 
employ a different technique that is taylored to explore the quantum-classical correspondence of the structure of 
eigenstates. It consists of representing the Hamiltonian (|I|) in new coordinates (w, v) in such a way that the effects of 
the boundary are transferred to an effective interaction potential between the two degrees of freedom (u, v). That is, 
the coordinates (u,v) are chosen so that the new wave function satisfies "flat" boundary conditions: \t(u, v) = at 
v = 0; d. For our rippled channel, this can be accomplished by the transformation 

u = x 

yd y , n 

V= d + a£(x) = 1 + U 

where e = a/d is a measure of the perturbation due to the ripple Jl7| ]. The Schrodinger equation in the new coordinates 
can be obtained from the covariant expression for a particle moving (in the absence of potentials) in a Riemannian 
curved space [ pj| , 

h 2 . , h 2 



A cov V(u, v) = £-g- 1 ' a d a g a Pg 1 ' 2 d p *(u, v). (3) 
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Here A cov is the Laplace-Beltrami operator, g is the metric and g a/3 is the metric tensor. Even though (||) is still the 
kinetic energy, the resulting differential equation takes a much more complicated form than the ordinary Laplacian. 
This is the price we have to pay when we transfer the effect of the boundaries onto the operator [the explicit form 
of the Schrodinger equation in (u,v) coordinates for the rippled billiard is given in pj]. Moreover, the coordinate 
representation of the canonical momentum has now the form J19|, 



Pn 



-%h 



d a + ^d a ln(g) 



(4) 



Substitution of this expression into the Schrodinger equation (||) determines the quantum Hamiltonian in covariant 
form PI 



H 



1 



2m, 



■g-^P^g^Ppg- 1 / 4 , a, /? = «,», 



(5) 



Substituting now the explicit expressions for the metric tensor g a @ and metric g, 



1 -tvj-u 
1 + C S 

-e<„ l+e 2 u 2 gg 
1+et (l+£«) 2 



(6) 



g = Det{g aP ) = [1 + e£uf 



(7) 



we get 



H = ^—g- 11 ^ \Pu(1 + eOPu + P v l+ / { V T U P V 
2m e I 1 + e£ 

} 5 1/4 



(8) 



where t; u = d^/du. Note that the rippled boundary is fully incorporated into the Hamiltonian operator. Since the 
boundary conditions in the new coordinates correspond to those of a "flat" channel, ^(u, v) = at v = 0; d. For our 
purposes it is convenient now to separate the Hamiltonian (8) as 



H = H° + V(u,v,P u ,P v ), 



where 



with [c. f. (4) and (7)] 



2m. 



iPu+Pv) 



p., = -ih 



d u + ^d u ln{g) 



, Py iflOy , 



(9) 



(10) 



(11) 



and V stands for the rest of the terms. 

We remark that this representation allows us to treat the original model of one free particle in the rippled channel 
as a ID model of two interacting "particles" identified with the two degrees of freedom u and v. Here H° is the 
Hamiltonian of two non-interacting particles. The eigenf unctions of H define the unperturbed basis in which the 
eigenstates of the total Hamiltonian H (u, v, P u , P v ) may be expanded. Such a representation turns out to be convenient 
for the study of chaotic properties of the model since one can use the tools and concepts developed in the theory of 
interacting particles (see, for example, [p0[). 

Since the Hamiltonian (||) is periodic in u, the eigenstates are Bloch states. This allows us to write the solution 
of the Schrodinger equation in the form v E'_b(u, v) — exp(ifcu)$fe(w, v) with $fc(it + 2ir, v) — $*;(«, v). For an infinite 
periodic channel the Bloch wave vector k = k(E) takes a continuous range of values, which we take to lie on the first 
Brillouin zone (— \ < k < |). 

By expanding <&&(«, v) in a double Fourier series the a-th eigenstate of energy E a (k) can be written as 
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E CZ m m k mn (u,v), (12) 

m— 1 n— — oo 

where 

<fimn( u > v ) =<u,v\m,n > k = 

(13) 

7f !/2g i/4 sin(m7ri;/d) exp [i(/c + n)u] 
are the eigenstates of the unperturbed Hamiltonian H (u, v) with corresponding eigenvalues 



4t(*0 = 



2m P 



{k + nf 



717T\ - 



(14) 



The factor tt -1 / 2 ^ 1 / 4 in Eq. (|l3|) appears from the orthonormality condition in the curvilinear coordinates [u, v). 
Thus, the exact eigenstates are expanded in a complete orthonormal basis satisfying the boundary conditions of the 
problem, i.e., a Galerkin series expansion 

So far, there has been no need to specify the ripple function £(x), except that it is periodic: £(x + 2n) = £(x). For 
concreteness, we consider from now on the dependence = cos(a;). The ripple profile is given by y = d + acos(x) 
where x, y, d, and a are dimensionless quantities. The latter are defined from the expression Y = D + A cos(2tt X / B) 
scaled to the period B. Therefore, x — 2tt^, y — 2n^, d = 2njt, and a = 2tt^, where X, Y, D, A, and B are 
dimensional quantities. 

This approach to solving the Schrodinger equation was used to calculate the energy-band structure and Husimi 
distributions of the rippled channel Jl0| , and to obtain the energy level statistics under the conditions of full and mixed 
chaos j22j. Specifically, in Ref. j2i| the energy level spacing distribution was shown to be Poisson, Wigner-Dyson, or 
intermediate between these two, for regular, globally chaotic, or mixed classical motion, respectively, in agreement 
with the well known RMT conjecture [£3[. Most relevant to our work here is the fact that these distributions were 
found to be the same for all values of k within the Brillouin zone, except k ss for which the parity symmetry 
should be taken into account. Thus, without loss of generality we fix the value k — 0.1. We shall also use in all our 
calculations the following values for the ampliutude a and width d of the rippled channel: a/27r = 0.06 and d/2ir = 1.0. 
These values produce representative global chaotic motion in the classical limit (see below) . 

In order to study the structure of eigenstates of the total Hamiltonian H (u, v) one needs, first, to choose a way of 
ordering the unperturbed basis in which to represent the Hamiltonian matrix Hi ^{k) =< I \ H \ V >/.. Specifically, 
we have to assign an index /, labeling the basis state | / >fc=| m,n to each pair of indecies (m,n) (note that, 
although the energy spectra is independent of the assignment (m,n) — > I, the structure of the eigenstates is not). 
The size of the Hamiltonian matrix is determined by the maximum values of n and m: —N max < n < N max and 
1 < m < M max . 

A natural assignment is the following one. Let us fix the lowest value of n {—N max ) and sweep all values of m (1 < 
m < M max ). This gives 1 = 1,2, ...,M max . Then do the same for n = -N max + 1, which gives M max + 1 < I < 2M max , 
and so on, till finally we have 1 < I < L max , where L max = (2N max + l)M max defines the matrix size, L max 

x L max . 

This rule results in a block structure of the Hamiltonian matrix, with block size equal to M max . Fig. 2 shows the 
central part of a 4030 x 4030 matrix with (N max , M max ) = (32, 62). Here we can see a number of blocks of size 62 x 62 
corresponding to n,n' — 0, ±1,±2 ± 3. In this representation the matrix is clearly a band matrix. The finite size 
of the band is due to the short range coupling between different blocks, namely, the strength of the matrix elements 
decrease as e l"-"'l, e l™-«'±i| ; r e l«-™'± 2 l ( see details in 0). 

The above way of ordering the unperturbed basis is typical; it corresponds to the "channel representation" since 
the index m labels a specific transverse channel (or mode) for the propagation of the wave through the billiard, 
see Eq. (|l3|). However, for our purposes it is essential to use the "energy representation", according to which the 
unperturbed basis is ordered in increasing energy, Ef +1 (k) > Ei(k). This defines a new rule / — >■ l new = l new {n,m). 
Fig. 3 shows the difference between these two ways of ordering of the basis, giving the normalized unperturbed energy 
E°(k) = E^(k)/E*=(^) 2 {n + k) 2 + m 2 with E* = as a function of I. 

A non-monotonic dependence of the energy occurs in Fig. 3a because the first rule l(m, n) described above produces 
minimum values in the energy for m = 1. In contrast, in Fig. 3b a linear dependence of the unperturbed energy on 
the index l new is shown, apart from large values of l new where the finiteness of the matrix results in the distortion of 
the spectra. The numerical obtained linear dependence, E = E/E* — ^l ne w, can be derived analytically from Weyl's 
formula N(E) = 4ppS valid for 2D billiards with area S (also valid for periodic 2D billiards) |2J . 
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A crucial point is that the eigenstates of the total Hamiltonian in the "energy representation" have a very convenient 
form for the analysis. The advantage of the "energy representation" is clearly seen in Fig. 4 where an arbitrarily 
chosen eigenstate is given in the two representations. One can see that in one case the eigenstate has a kind of regular 
and extended structure, while in the other, the eigenstate is compresed around the unperturbed state whose energy is 
close to the energy of the perturbed state. In the latter case one may use a statistical approach to describe the global 
properties of such eigenstates, see Refs. |2(],|25|,^6|] . Specifically, we may characterize such eigenstate by introducing 
an envelope around which the components are expected to fluctuate in a pseudo-random way. We stress that by 
using this energy ordering one can relate the global form of eigenstates in the energy representation with its classical 
counterpart, see below. 



III. CLASSICAL REPRESENTATION OF THE MODEL 



Clearly, the type of motion of a particle in the ripple billiard depends on the values of the geometrical parameters, 
namely, on the ripple amplitude a and the average width d. If the channel is narrow, d <C 2ir, the Poincare section 
reveals a large resonance island surrounded by the typical Poincare-Birkhoff structure . This island is formed 

by the librational motion along the x direction in the neighborhood of a; = 0. On the contrary, for wide channels 
(d > 2tt) global chaos occurs even for small amplitudes due to a strong overlap of resonances. In this work we shall 
limit ourselves to the study of global chaos, i.e., wide channels. 

The condition of global chaos can be derived analytically for the case of small amplitude a <C d where the following 
approximate map is valid |B7fl, 



../,■■>-. 



= a„ 



■ 2a sm(x n ) 
2d cot(a„+i) 



(mod 2tt), 



(15) 



Here x n is the position of a particle and a n is the angle between the x-axis and the velocity of the particle, right after 
the n— th collision with the upper wall. The standard linearization around fixed points of period one yields, 



5a n+ \ = Sa„ + k sin(Ax„) 

Ax n +i = Ax n + 5a n +i, (mod 2tt), K = 2da/n 



(16) 



where the angle 5a and the position Ax, measured from the fixed point, play the role of the action and angle variables, 
respectively. Eq. ( |l6| ) has the form of the standard map with k as the nonlinear parameter. Hence, Chirikov's 
overlapping criteria predicts the onset of global chaos for k ~ 1. This was confirmed by computing numerically the 
actual path of the particle as it travels along the channel. Specifically, for k > 1 the Poincare section (a n ,x n ) shows 
global chaos for wide channels. In this paper we shall consider a wide channel with small ripple amplitude, specified 
by the parameters (a/2ir, d/2ir) = (0.06, 1.0), therefore, k w 1.5. 

Let us now write, as in the quantum description, the classical Hamiltonian in the (u, v) coordinates so that the 
effects of the rippled boundary can be incorporated as a kind of coupling between two degrees of freedom. The 
classical Hamiltonian is 



H = 



1 



2m, 



(17) 



which can be obtained from the quantum Hamiltonian (j8|) by commuting all momenta and coordinates |2q] . 

For completeness, the canonical transformation between (x, y, P x , P y ) and (u, v, P u , P v ) are given by Eq. (0) together 
with P u = P x + yP v /(l + e£), and P v = (1 + e£)Py Inserting the metric tensor Eq. (||) into ( |l7|) and regrouping 
terms, the Hamiltonian can be written as H = H + V, where 



and 



V 



2m e 



1 



2m, 



(Pu+Pv) 



2Z + e(e 



1+6^ 



PuPv 



(18) 



(19) 



V represents the interaction between the two "particles" u and v (here £ u = d^/du, and e = a/d). 
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IV. STRUCTURE OF EIGENSTATES AND LDOS: DEFINITIONS. 



Once the matrix Hi ii(k) has been diagonalized, its eigenstates ^ a (k) = are a ^ so re-ordered in energy: 

E a+1 > E a . We adopt the convention that the greek superindex (latin subindex) denotes the exact (unperturbed) 
state. The amplitudes Cf(k) form the "state vector matrix". The elements along the row a of this matrix are the 
components of the a th eigenstate in the representation of the unperturbed energy-ordered basis. Correspondingly, 
the elements along the column I give the unperturbed state I expanded in the energy ordered perturbed basis. For 
our purposes, we examine the matrix wf = |C"(fc)| 2 (see Fig. 5) which plays the central role in our approach. 

The rows of the matrix wf show how a specific eigenfunction (eigenstate \a >) is expanded in the unperturbed 
basis \l >. Since the average "band width" of wf smoothly depends on the index a (equivalently on the energy E a ), 
it is convenient to average wf over a small range of eigenstates. In this way, one can obtain the so-called structure of 
eigenstates (SES) in the unperturbed basis. In what follows, we are interested in the SES in the energy representation 
which is defined as p9[, 

W(E? \E a )=J2 < S(E$ - E a '), (20) 

a' 

where wf = wf /N with TV as the number of eigenstates \a' > in the vicinity of a given a. The structure of eigenstates 
W(E? | E a ) gives the dependence on the unperturbed energy Ef for eigenstates with total energy close to E a . This 
sum can be approximated by an integral that involves the density of states p{E a ) providing the transition from the 
basis to the energy representation. This function SES plays the key role in the determination of important physical 
quantities, such as the distribution of occupation numbers in chaotic closed systems (see, for example p|, p0| , |30"| ) . 

Analogously, one can analyze the columns of the matrix wf which give information about a given unperturbed state 
\l > is spanned over exact eigenstates \a >. Therefore, one can define the following quantity 

co(E a \Ef) = Y,<>KE a -Ef), (21) 
v 

which is well known in nuclear physics as the strength function, and in solid state physics, as the local density of states 
(LDOS). In (pl| ) the summation (average) is done over N basis states \V > around a fixed state |Z >, wf = wf, /N. 
Therefore, this quantity is considered as a function of the energy E a . It shows how a given unperturbed state | <f>i > 
is coupled to the exact states due to the interaction V. The width of this function (spreading width) determines 
the energy range associated with the "decay" of a given unperturbed state into other states when the interaction is 
switched on. 



V. CLASSICAL ANALOG OF THE STRUCTURE OF EIGENFUNCTIONS 

An essential point is that both, the structure of eigenstates and the local density of states, have well defined 
classical analogs p9[ . Let us first start with the SES. Since Cf —< i S a \ <pi > as a function of I is the projection of 
the perturbed state onto the states of the unperturbed system, the classical counterpart of wf =| Cf | 2 as a function 
of energy Ei can be defined as the projection of the total Hamiltonian H onto the unperturbed one H , where 
H = H° + V with V the perturbation |29|. This projection can be numerically done by substituting the trajectories 
<i>(i) = (x(t),y(i),p x (i),Py(t)) generated by the Hamiltonian H with energy E onto the unperturbed Hamiltonian 
H°. Since the unperturebed energy E (t) along these trajectories varies in time, it fills the so-called energy shell 
characterized by its width AE. Therefore, for chaotic total Hamiltonians H, the classical analog W c i(E° | E) of the 
quantum SES can be easily obtained from E°(t). 

In Fig. 6a we show the energy of the unperturbed Hamiltonian H° as a function of time after the substitution 
of a single chaotic trajectory <f>(t) generated by the perturbed Hamiltonian H, into H°. The classical distribution 
W c i(E° | E) (see Fig. 6b) is constructed from E°(t), averaged over a sufficiently long time. Alternatively, since the 
dynamics of H is chaotic, the same distribution can be obtained by averaging over many different orbits at a shorter 
time. For concreteness, in what follows we chose E = 1. 

It is important to understand the origin of the shape of the classical distribution, see Fig. 6b, in view of its 
correspondence to the quantum SES. For this, let us examine in detail characteristic time intervals of the plot of E (t) 
together with plots of P 2 , P 2 , and v, see Figs. 7a-d. From Fig. 6a it is clear that the main contribution to the central 
peak of the distribution W c i(E° | E) comes from the continuous time intervals where E° w 1. These intervals in turn 
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correspond to P 2 w 2 (the maximum value) and to P 2 ~ (see Figs. 7b-c). During these time intervals the particle 
travels almost parallel to the x axis (u axis). In fact, from the plot of v(t) (Fig. 7d) one can see that the particle hits 
only once the bottom boundary during a relatively long time. The contribution to the central peak can be also be 
seen directly from the expressions ( |l8| ) and ([l9|). Indeed, since P 2 w 0, we have E°(t) w ^P 2 ~ E = 1. 

In the same way one can understand the origin of two side-band peaks in the form of W c i (E° | E) , clearly seen in Fig. 
6b. The data reported in Fig. 8a-d show that for small time intervals e. g., (0.212 < t < 0.216 and 0.231 < t < 0.234) 
the value of E° is practically constant (E° w 1.12 and E° « 0.88 resepectively) . These values of E° correspond to 
the right and left peaks in Fig. 6b. Contrary to the data in Fig. 6 explaining the origin of the central peak, in the 
two time regions of Fig. 8 the value of P 2 is nearly zero and P 2 is nearly maximum (see Fig. 8b-c). In these regions 
the motion of the particle is almost perpendicular to the x-axis (it-axis). 

The form of the profile £(it) = cos(m) in our ripple channel results in two periodic orbits of period one: the stable one 
for P u — 0, u = 0, and the unstable one for P u = 0, u — tt. Putting these values for P u and u into expressions (18-19), 
we find the values for P 2 . This allows us to obtain the value of H° for these two periodic orbit, E° = E(\ + e) 2 w 1.12 
and E° = E(l — e) 2 w 0.88, respectively for u — and u = tt. Thus, the left (right) peak of the distribution 
W c i(E° | E) is formed by trajectories dwelling near the unstable (stable) periodic orbits of period one. 

VI. CLASSICAL ANALOG OF THE LDOS 

In analogy with the quantum LDOS, the classical LDOS distribution u> c i(E | E°) is constructed by projecting 
the dynamics generated by H° = E° onto the Hamiltonian H = H° + V. Here, because H° is integrable, several 
(regular) trajectories of H have to be substituted into H and then averaged in time, in contrast with the classical 
SES where only one (chaotic) trajectory over a sufficiently long time is needed to form the distribution. The result of 
this procedure is exemplified in Fig. 9a where E(t) is shown for 36 regular trajectories, whose initial conditions were 
chosen to lie on a mesh of points distributed uniformily in the plane (P u , u) (v is fixed and P v is determined by the 
constant energy E° = 1). The classical LDOS distribution us c i(E | E°), see Fig. 9b, was obtained from E(t) (Fig. 
9a). 

In the construction of Fig. 9a, 6 initial conditions for u (taken from the interal [— tt,tt]) were used for each of the 6 
initial values of P u [from the interval (— V2E°,V2E t) )]. Hence the appearance of 6 similar sets formed by 6 small time 
intervals. The very first time interval ([0,0.25]) yields E = 1 and is produced by a trajectories with initial condition 
P u ~ — v/2 and u — —tt, i.e., a grazing trajectory. The next interval [0.25, 0.5], a bar code-like structure extending over 
the whole energy shell, is produced by the initial condition P u = — 0.6V2, u — —tt, and the third interval [0.5,0.75], 
the noise-like one, is produced by the initial condition P u = —0.2^/2, u = — tt. The next three intervals are produced 
by mirror initial conditions (P u = 0.2-\/2, 0.6-\/2 , and ~ \[2 with u — —tt). The initial conditions for the next sets 
are the same as for the first set, except u is shifted subsequently by 2tt/5. Any trajectory reveals one of these three 
types of behavior of E(t); E = 1, bar code, and noise-like structure. The averaging over these time intervals yields 
the distribution shown in Fig. 9b, quite different from the classical SES (Fig. 6) although it also has the central peak 
and two other peaks on each side. 

VII. QUANTUM-CLASSICAL CORRESPONDENCE: GLOBAL PROPERTIES 

Having understood the dynamical origin of the main features of the classical distributions W c i(E° | E) (SES) and 
uj c i(E I E°) (LDOS), we now compare them with the corresponding quantum quantities. Although the SES for a 
given a, see Eq. (20), is actually the average over a set of perturbed states in the neighborhood of a, it is convenient, 
first, to examine typical individual states in different energy regions. In order to compare the quantum and classical 
behavior, we use the quantum energy of interest E q (where E q is E a in the case of SES and in the case of 

LDOS) to construct the classical distributions (E c i = E q ). The ratio of the de Broglie wavelength A = 2tt^J 2m~E; *° 
some characteristic length of the billiard may be used as a semiclassical parameter (recall that there are two length 
parameters defining the flat channel, D and B; in this work we have used D = B). For example A/D = 2Trti e ff << 1, 
where h e f f = E , tells us how large the energy of the eigenstate should be for its wavelength to be many times 

smaller than the width of the channel. Substitution of Weyl's formula N(E) = B 2 ^l" E into the above expression for 
A gives a more useful formula: 
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A/£> = /fj = 2vK eff . (22) 

Here we only need to know the number of the eigenstate (as oppossed to the energy) to determine the ratio A/D. 

Fig. 10 shows wf for four representative states (a = 111, 712, 1362, and2002) of different energy regions as a function 
of the unperturbed energy scaled to E a . The classical SES is also shown. As it may be expected, these figures 
show that as S, e // becomes smaller the number of unperturbed states needed to form the pertubed state increases. 
Specifically, level states around a = 111 are still far from the semiclassical regime (A/D ss 0.17); the details of the 
rippled boundary cannot be resolved and, thus, there is only a weak mixing of the unpertubed states. In contrast, for 
the state a — 2002 its de Broglie wavelength is a smaller fraction of the width of the channel (A/D ss 0.04) and the 
number of participating components is larger. 

It is important to note that although the number of unperturbed states needed to construct the perturbed state 
increases as the energy increases, they all fall within the range of energies determined by the classical distribution. 
This range defines the energy shell of the eigenstate under consideration. On the other hand, it is clear that even 
the level a = 2002 does not seem to have much of a resemblance with the classical distribution. For this and similar 
states, we have calculated the distribution of fluctuations arround the classical curve and found that they deviate 
strongly from a Gaussian distribution as predicted by RMT considerations. This may be related to the fact that the 
de Broglie wavelength of level 2002 is still not much smaller than the amplitude A of the ripple (A/ A ss 0.65, recall 
A/D=0.06). That is, even though A/D may be small, a more relevant parameter to characterize the semiclassical 
regime for particular eigenstates is A /A. Nevertheless the strong fluctuations exemplified by the level a = 2002 (Fig. 
lOd) can be smoothed out by averaging over a range of perturbed states. This actually corresponds to the definition 
@ of the structure of eigenstates W(E? \ E a ). 

Remarkably, as Fig. 11 shows, such an average does reveal a good quantum-classical correspondence even though the 
fluctuations are strong. Specifically, the global shape of the quantum SES displays a three-peak structure, much like 
the classical one. In addition, the energy spread of the quantum SES agrees very well with the classically determined 
energy shell, even in its slight assymetry (0.7 < E° < 1.45 with the center at E° = 1). 

Similarly, in Fig. 12 we compare the quantum and classical LDOS. The quantum-classical correspondence appears 
to be even better than for the SES in the sense that the fluctuations are smaller. 

VIII. LOCALIZATION AND NON-ERGODICITY 

To understand the origin and importance of the strong fluctuations of the SES around the classical counterpart (see 
Fig. 11), we analyze the structure of individual eigenf unctions. In Fig. 13 three consecutive typical eigenfunctions 
(a = 1952, 1953, and 1954) are shown. The difference between the state a = 1954 and the other two is clearly 
qualitative. More specifically, while the states a = 1952 and 1953 are extended (in energy) states, constituted 
by practically all unperturbed eigenstates within the energy shell, the state a — 1954 is mostly unperturbed; it 
is extremely localized in the energy shell. By neglecting all small amplitude components surrounding the main 
component (see Fig. 13c), we can determine the unperturbed state, defined by a pair (m,n), that most closely 
resembles the perturbed state. We find that this always corresponds to the lowest values of the transversal mode 
to. This fact can be understood by the following physical argument. Consider an eigenstate of the flat channel 
$ nn (X : Y) cx sm(mTrY/D) exp(iKX) with energy E° = ^(K% + K%), where K v = ^ = f^. Turning on the 
perturbation (flat to rippled channel) will affect the high energy unperturbed states differently depending mainly on 
the value of A y . For example, for to = 1, the ratio A y /A is 2D /A ps 33, which is so large that the state cannot 
"see" the ripple and thus will remain essentially unperturbed. In contrast, for unperturbed states with the same (or 
about the same) energy but with large values of to (say, m = 62 = M max and correspondingly small k x ) their A y 
is sufficiently small, compared to the amplitude of the ripple (A y /A « 0.5), that the rippled boundary produces a 
strong mixing of unperturbed levels. The resulting perturbed state will consist of many components, extended over 
the energy shell. By the same token, it is expected that unperturbed states with intermediate values of to will turn, 
after turning on the perturbation, into some intermediate states, the so-called sparse states [32], made up of fewer 
components over the energy shell. 

The existence of the extremely localized (in energy) states manifests itself in the structure of the Hamiltonian 
matrix Hi ji in the channel representation (Fig. 2) to be discussed below. It is easy to find that such wavefunctions 
differ only slightly from the plane waves c\P m n (x, y) with small m, proper of the flat channel. In Fig. 14a one localized 
(in energy) state is shown in configuration representation. In contrast, an extended (in energy) state is presented in 
Fig. 14b. 
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Eigenfunctions like that of Fig. 13b for a = 1953, are somewhat intermediate between the extended and localized 
states. 

In order to characterize quantitatively the eigenfunctions, we compute various localization measures. The first one 
is the so-called entropy localization length Ih, 



l H = cxp [- (H - Ugoe)] ~ 2.08 exp {-H). 
Here Ti. stands for the Shannon entropy of an eigenstate in a given basis, 



n 



N 

E 

i=i 



wf Inw?, 



(23) 



(24) 



and TLgos is the entropy of a completely chaotic state which is characterized by gaussian fluctuations (for N — > oo) 
of all components Cf with the same variance < wf >— l/N. The latter property occurs for completely random 
matrices belonging to a Gaussian Orthogonal Ensemble (GOE). Defined in this way, the quantity Ih gives a measure 
of the effective number of components in an eigenstate. For example, the eigenstates of Fig. 13 have, respectively, 
l n = 745.9, 232.9, and 2.7. 

The second quantity, U pr , which gives another measure of the effective number of components in an eigenstate, is 
expressed via the inverse participation ratio V, 



I , 



V 



3 

V' 



with 



V 



N 

E 

1=1 



(25) 



(26) 



where Vgos ~ 3 is chosen in order to get li pr = N in the GOE limit case. 

The above two definitions of "localization lengths" are the most frequently used when describing global structure 
of eigenfunctions. One should note that these quantities give an effective number of large components, independently 
on the location of these components in the unperturbed basis. 

Additional information about the structure of eigenfunctions can be obtained from the "width" or mean square 
root l a of an eigenstate, computed as 



n 



$>? [l-n c (a)Y 



1=1 



1/2 



(27) 



where n c = I wf determines the centroid of an eigenstate in the unperturbed basis. 

To get a complete panorama, in Fig. 15a-c we plot these three measures as a function of a corresponding to the 
eigenstate \a >. The strong fluctuations of all localization measures are evident in these figures. We can see that 
neighboring high levels may have drastically different localization measures, in agreement with the discussion above 
about the existence of the three types of states: localized, extended, and sparse. These figures gives us information 
about the relative number of each type to be found in a given energy range. 

Comparison of the width l a with Ih and U pr gives the possibility to detect the so-called sparsity of eigenstates. 
Indeed, small values of the ratio In/la (or h pr /la) indicates that there are many "holes" in the structure of eigenstates, 
therefore, such eigenstates are sparse p2[ . A detailed analysis shows a dominance of sparsed eigenstates. As for the 
centroid n c — a small fluctuations (observed in Fig. 15c) indicate that the interaction strength V is relatively weak, 
compared to the unperturbed part Hq. 

The data of Fig. 15 show the existence of some pattern for all localization measures, visible as clusters of points 
in the lower part of the plots. These patterns are much more pronounced for the same quantities computed for the 
LDOS, see Fig. 16. In order to understand this unexpected phenomenon, let us take a closer look at l a in the range 
1450 < I < 2100 and examine the structure shown in Fig. 17 in detail. Specifically, we choose some eigenstates from 
three branches of points (see marked states in Fig. 17). From each of these three sets, we take the 1 st , 7 , 13 th , 
and 19 th states, counted from the bottom of the figure, and plot them in Fig. 18. Inspection of Figures 18a-18d 
clearly demonstrates that there is a kind of regularity in the structure of the LDOS: the same type of states appear 
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repeatedly, almost periodically as a function of the basis number I. These figures show the repetition of extremely 
localized states (Fig. 18a) and of states with different degrees of sparsity (Figs. 18b-18d). The repetition of extended 
states, corresponding to larger vaues of l a are not shown for economy of space but their existence is clear, as can 
be inferred by extrapolating the branches of Fig. 17 to higher values of l a . The physical origin of all these types 
of states (localized, sparse and extended) was explained above, and their appearance can be decoded by examining 
the structure of the "channel representation" of the Hamiltonian matrix (Fig. 2). A detailed inspection of this matrix 
(see also Eq. (26) of Ref. p0||) shows that the coupling between unperturbed states depends strongly on the values 
of the index to, labeling the transversal modes of the flat billiard. An unperturbed state specified by a large value of 
to (an to close to M max — 62) couples strongly to several other unperturbed states. In contrast, a state with to = 1 
has practically no coupling to other states. In particular, the extremely localized states, corresponding to the first 
position on the left line of each brach of Fig. 17 occur because of the negligible coupling of the diagonal elements of 
the Hiji matrix with to = 1, the states on the second position of each branch occur for m = 3, and so on with to odd. 
Similarly, states on the right side of the branches result from elements of the Hiji matrix with even values of to. 

This structure is expected to prevail at all energy ranges since in any sufficiently large range of energies there are 
unperturbed states with all values of m in [1, M max \. Even deep in the semiclassical regime, extremely localized and 
sparse states will appear but less and less frequently since the energy differences between states of the same type 
increases with energy. Consequently, the strong fluctuations appearing in the SES and LDOS distributions (Figs. 11 
and 12) will tend to dissapear as h e ff — ► 0. 

In the definitions of SES and LDOS, the averages were performed over the appropriate range of energy in order to 
take into account once all the various types of states (Figs. 10 and 11 were calculated this way). 

We can understand the origin of this type of regular structure by analyzing the Hamiltonian matrix H^i' , see Fig. 
2. A detailed inspection shows that for the unperturbed states (diagonal matrix elements), the coupling with other 
unperturbed states depends on the value of the index to which labels the number of transversal modes in the billiard 
[see Eqs. (11-12)]. An unperturbed state specified by a large m (an m close to M max = 62, for any n) shows a 
strong coupling to several other unperturbed states. In contrast, a state with to = 1 (for any n) shows practically no 
coupling to other states. 

We can see then that the extremely localized states, corresponding to the first position on the left hand side of 
each set of points in Fig. 16 occur because of the negligible coupling of the diagonal elements of the Hij> matrix with 
to = 1. Similarly, the states on the second position (on the left hand side) occur when to = 3, and so on with odd 
to. States on the right hand side of each set of points result from elements of the Hi^i matrix with even values of 
to. This structure prevails for all energies due to the existence of unperturbed to = 1 states. However, the intervals 
between the occurence of to = 1 unperturbed states increase as the energy increases. Even as the semiclassical limit is 
approached, extremely localized states will appear, except less and less frequent. Consequently the strong fluctuations 
appearing in the SES and LDOS distributions (Figs. 10 and 11) will tend to dissapear. 

IX. CONCLUDING REMARKS 

The quantum-classical correspondence for the chaotic motion of an electron in a periodic billiard was analyzed in 
terms of the classical analogues of the structure of eigenstates (SES) and local density of states (LDOS). 

To construct the classical counterparts of the LDOS and SES, we first changed (via a canonical transformation) 
to some new curvilinear coordinates, where the new Hamiltonian incorporated the effects of the boundary into an 
effective interaction potential. Then the original system of a free particle colliding within a 2D rippled channel becomes 
equivalent to a ID model of two interacting "particles", identified with the new coordinates. This procedure allowed 
us to study the chaotic properties of this system, using tools developed recently to analyze the role that interactions 
between particles play in the onset of quantum chaos. This example is quite instructive since it shows the relevance 
of "one-body chaos", which is due to an external potential (or boundary conditions), to "two-body chaos", which 
results entirely from the interaction between particles. So far, these two mechanisms for producing chaos have been 
treated completely independently 

Using this approach, we found that in the case of strong chaos, the classical analogs of SES and LDOS agreed 
remarkably well with the global shapes of the quantum quantities. This correspondence enabled us to explain the 
main features of the SES and LDOS in terms of the underlying classical motion of the system. Specifically, we 
found that classical unstable and stable periodic orbits of period-one give rise to the two pronounced side-band peaks 
observed in the SES and LDOS. 

On the other hand, we found that the statistical properties of eigenstates and individual LDOS differ qualitatively 
from those prescribed by the standard random matrix theory. Namely, in the case of strong chaos and deep in the 
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semiclassical region (high energy), one expects the components of individual eigenstates to fluctuate in a statistically 
independent way around the mean, the envelope of the SES in the unperturbed energy basis. This expectation was 
confirmed when studying the structure of eigenstates of complex atoms and nuclei in the mean-field basis p4] , |35|] . 
In contrast, we have found that in our present model the deviations turn out to be extremely strong and not fully 
statistically independent. 

In connection with this, a detailed analysis revealed quite unexpected regularities in localization measures (such as 
the inverse participation ratio) characterizing the eigenstates and individual LDOS. We remark that these regularities 
do not disappear as energy increases. In particular, even for high energies, one can find eigenstates that are strongly 
localized in the unperturbed energy representation (i.e., slightly perturbed plane waves). This occurs because energy 
is not the only semiclassical parameter in a 2D electron waveguide. Clearly, two unperturbed states of similar (or 
equal) energy values but different transverse mode numbers will react differently to a small perturbation. That is, in 
contrast to high-mode unperturbed states, the transverse wavelength A y for low-mode states can be many times larger 
than the amplitude of the perturbation and hence will remain essentially unperturbed. 

Thus, for any value of energy (equivalently, for any small value of K), there are non-ergodic states. It should be 
stressed that here we are discussing the ergodicity in the energy shell, which is determined by the classical motion 
of a particle in terms of the classical SES. Since the width of the energy shell is always finite due to the finite range 
of the interaction, the ergodicity of an eigenstate means that this eigenstate fills the whole energy shell with random 
(Gaussian) fluctuations of its components around the smooth energy dependence defined by the classical SES. 

Naturally, with an increase of the energy, the relative number of localized (in the energy shell) eigenstates tends to 
zero. In this sense, there is no contradiction with the onset of quantum ergodicity in the classical limit. However, this 
limit turns out to be achieved very slowly. Therefore, we remark that from a physical point of view, it is important to 
further the study of the statistical properties of eigenstates in the deep semiclassical regime, which may be different 
from those in the strictly classical limit {% = 0, a mathematical concept). 
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FIGURE CAPTIONS 



FIG. 1 Geometry of a rippled billiard. 

FIG. 2 Central part of the 4030 x 4030 Hamiltonian matrix H hV {k): N max = 32, and M max = 62. The 62 x 62 
blocks corresponding to n,n' — 0, ±1,±2, ±3 are shown. 

FIG. 3 (a) Unperturbed energy spectra as a function of I. (b) Unperturbed energy spectra ordered in energy as a 
function of l new . 

FIG. 4 Example of eigcnfunctions for different ways of odering the unperturbed basis, (a) the state a = 2122 as a 
function of I; (b) the same state as a function of l new . 

FIG. 5 Left upper part of the matrix wf corresponding to N max = 32, and M max = 62. 

FIG. 6 a) Energy of the unperturbed Hamiltonian E°(t) as a function of time (in arbitrary units) for E = 1. Note 
the intermittent structure (from a noise- like to a stable-like) in E°(t). b) Classical distribution W c i(E° \ E) 
constructed from E°(t). 

FIG. 7 Energy E°, and P%, P%, v as a function of time (in arbitrary units) on a short time scale, after the substitution 
of a trajectory <j>(t) £ H onto H° for E = 1. 

FIG. 8 The same as in Fig. 7 for a different time interval. 

FIG. 9 a) Energy E(t) of the full Hamiltonian H as a function of time (in arbitrary units) after the substitution of 
36 trajectories <j>°(i) £ H° onto H for E° — 1. b) Corresponding classical LDOS ui c i(E | E°) constructed from 
E(t). Note that each trajectory contributes in a different way to the classical distribution as can be seen in Fig. 
9a, for details see the text. 

FIG. 10 Individual eigcnfunctions wf in the energy representation and the classical distribution W c i(E° \ E) (thick 
curve) as a function of the energy Ef jE a for a= a) 111, b) 712, c) 1362, and c) 2002. The effective Planck's 
constant h e ff is a)0.027, b) 0.011, c) 0.008, and d) 0.006, respectively. The entropy localization length (defined 
in Sect. VIII) l H is a) 23.6, b) 158.9, c) 483.6, and d) 709. 

FIG. 11 Structure of eigenstates (SES) in the energy representation W(Ef \ E a ) and its classical counterpart 
W c i (E° | E) (thick curve). The average for the SES is taken over the interval 1880 < a < 2025. 

FIG. 12 LDOS uj(E a \ Ef) and its classical counterpart uj c1 {E \ E°) (thick curve). For the LDOS the average is 
taken over the interval 1900 < Z < 2000. 

FIG. 13 Typical high energy eigcnfunctions for a) a=1952, b) a=1953, and c) a=1954. The localization measures 
(defined below) are: a) Z ff =745.9, Z ipr =529.7; b) Z ff =232.9, Z tpr =45.4; and c) Z#=2.7, h pr =3.2. 

FIG. 14 Left: a localized and an extended eigenfucntions in the basis representation for (a) a = 2082 (extremely 
localized state with corresponding m — 2), (b) a = 2083. Right: the same eigenfunctions in configuration 
representations. 

FIG. 15 Localization measures for eigenfunctions in dependence on a for exact a th eigenstates. (a) Entropy local- 
ization length (b) l ipr defined through the inverse participation ratio, (c) the mean square root l a , and (d) 
the centroid n c . 

FIG. 16 Localization measures for the LDOS with I standing for the unperturbed state \l >. Same localization 
measures as in Fig. 15 are presented: (a) lu, (b) hp r , (c) l a , and (d) n c . 

FIG. 17 Detail of Fig. 16c for the range 1450 < I < 2100. The numbers indicate specific states which are shown in 
Fig. 18. The states are chosen from 3 different sets of states, specified by almost the same values of l a . 

FIG. 18 Structure of individual LDOS marked in Fig. 17. One can see the similarity in the structure of the LDOS 
taken from different sets, but with the same values of l a . 
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